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It is shown that nuclear level densities in a finite space are described by a continuous binomial 
function, determined by the first three moments of the Hamiltonian, and the dimensionality of the 
underlying vector space. Experimental values for ^^Mn, ^^Fe, and ®°Ni are very well reproduced by 
the binomial form, which turns out to be almost perfectly approximated by Bethe's formula with 
backshift. A proof is given that binomial densities reproduce the low moments of Hamiltonians of 
any rank: A strong form of the famous central limit result of Mon and French. Conditions under 
which the proof may be extended to the full spectrum are examined. 
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The vast majority of studies on level densities rely on 
Bethe's formula M, 



PB{E,a,A) = 



/27r 



^y/ia{E+A) 



12 (4a)i/4(^ + A)5/4 



(1) 



Although only two parameters are involved for a given 
nucleus — and a third, the spin cutt-off, for work at fixed 
angular momentum J — their dependence on mass num- 
ber, excitation energy and shell effects makes their de- 
termination delicate (see |g,g] for recent work and earlier 
references) . The problem is that the validity — and hence 
the success — of the formula goes well beyond the inde- 
pendent particle assumption made in deriving it, and one 
is left wondering about the true meaning of the parame- 
ters. The Shell Model Monte Carlo method provides an 
alternative approach to level densities [Q^^ whose reli- 
ability is now established Q. The problem is that the 
calculations are hard. 

To combine the simplicity of Eq. (|l|) with a 
parametrization of clear microscopic origin, we propose 
a continuous binomial form to describe the shell model 
level densities. It will be shown to do well with experi- 
mental data, to reproduce strikingly Eq. (1^) over a wide 
range of energy, and to provide a strong form of the cen- 
tral limit theorem (CLT), generalizing the famous result 
of Mon and French B (MF) for Hamiltonians of arbitrary 
rank. 

Consider a system of m particles moving in D orbits, 
spanning a space of dimensionality d. To specify a bino- 
mial density pb(x,N,p,S), three parameters are needed: 
N, the effective number of particles, the asymmetry p, 
and an energy scale e. The span (distance between low- 
est and highest eigenstates) , centroid E'c, variance a"^ and 
the adimensional energy variable x are given by 

S = Ne, Ec = Npe, a^ = Npqe^, x^^, (2) 
Where p + q = 1. Calling x ^ I — x, the density is 



Pb{x,N,p,S) 



: p^^q^^d 



T{N + 1) 



N 



T{xN + l)r{xN+l) S' 



(3) 



which reduces to a discrete binomial, ( ) if x = n/N 
ns/S, with integer n. 

It is often convenient to introduce A = p/q. Then 



p^^q^^d ■ 



\xN 



(1 + A) 



N 



X^"do 



do{i + xy\ (4) 



where do is the number of states at a; = 0. To deter- 
mine N, p and S (or e), Eqs. (0) for S and Ec cannot be 
used , since the spectrum is not known. We have to rely 
instead on the moments of the Hamiltonian 7i, i. e., av- 
erages given by the traces of H^ , to be equated with the 
corresponding moments oi pb{x, N,p, S), which for low K 
are the same as those of the discrete binomial (Eqs. (|ll|)). 
The necessary definitions and equalities follow. 

d-itr(H^) = in''), E, = in'), Mk = an - E,f) 



M^, Mk = 



M 



K 



K 



72 = Mi - 3 = 



a 
6pq 



Npq 



71 = Ms == 
d 



q~p 



do{l+p/qf. (5) 



N and p can be extracted either through the equations 
for 7i and 72 , or through those for 71 and d. The former 
option is unambiguous, and it has the advantage of warn- 
ing us that the binomial form is doomed if 72/7? > 1- 
The latter (which we adopt here) is simpler, and phys- 
ically cogent for the natural choice do=l, which locates 
the ground state at a; = 0. Once N and p are known, 
5* = [Na^/pqy/-^ follows. The centroid Ec provides the 
energy reference. With the simple choice, the predicted 
ground state is at Eq = —Sp, which may not coincide 
with the exact value, usually taken as origin. Therefore 
a shift A may be necessary, as in Eq. ([I]). 

Let us apply these prescriptions to ^^Mn, ^^Fe and 
^"Ni, for which data are available M. There is also a 
Monte Carlo calculation that does very well for the first 
of these nuclei m . The space chosen by Nakada and Al- 
hassid is the pf + gg/2 shells. The dimensionalities are 



given by d = {^){^) where D = 30 and z, n are the 
numbers of active protons and neutrons (e. g. z = 8, 
n — 12 for ^°Ni). In general, to proceed, we would 
have to calculate the three lowest moments of the Hamil- 
tonian, a feasible task even in enormous spaces. Here 
we simplify matters by making first the neutral choice 
p = 0.5, TV = lnd/ln2 = (42,44,49) for the three nuclei 
in the order above. Then, we bypass the calculation of 
a^ and estimate S directly . A first approximation comes 
from the single particle energies at (0, 6, 7, 7, 10) MeV 
for (/7/2,P3/2,/5/2,Pi/2,59/2) respectively, for which the 
differences between highest and lowest states come at 
(138,148,170) MeV for (^5Mn,56Fe,60Ni). Correlations 
increase these numbers by an amount that ranges from 
10 MeV for ^^^Ni to 20 MeV for ^^Cr, according to ex- 
act calculations in the pf shell |Q . Hence the estimate 
5'=(153 ± 5, 133 ± 5, 185 ± 5) for the trio. Fixing N and 
allowing S to vary within the error bars, the remaining 
uncertainty comes from the position of the ground state, 
which is also allowed to vary within zbe/2 « 1.5 MeV. 
Numbers that fall very confortably within the assigned 
error bars: S = (150,165,185), A = (1.15,-0.40,0.20) 
MeV, define binomial densities that agree very well the 
experimental ones ||ll| ( where the error bars are appar- 
ently smaller than the original ones in 0): 
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FIG. 1. Experimental, binomial (Eq. (|3|)) and Bethe 
(Eq. (0)) logarithmic level densities in ^^Mn and ^"^Fe. For 
parameters see text. 



Fig. w gives the experimental points for the first two 
nuclei, the binomials (whose parameters have been es- 
timated above), and fits to the binomials using Bethe's 
formula yielding (4a, A) = (21.5,-1.5) and (21.3,-3.4) for 
^^Mn and ^^Fe respectively. 

The near identity between binomial and Bethe forms 
extends to ~ 40 Mev. Therefore, the inmense experience 
accumulated by fitting Eq. (mj to the data can be reana- 
lyzed in terms of Eq. (ph whose phenomenological poten- 
tial has been illustrated by the simple exercise above. Its 
efficiency in dealing with rigorously calculated moments 
will be demonstrated in |13|. 



The near identity between Bethe and binomial forms 
has also a mathematical significance that is discussed in 
the paragraph that follows the one containing Eq. ( p^ . 

For ^°Ni in Fig. 0, instead of showing Inp^ — which 
fits obviously as well as in Fig. [I] — a different check is 
proposed. As our finite space can provide reliable den- 
sities only in a finite interval, to have an idea of its 
range, the space has been enlarged by adding the ^3/2 
and si/2 shells, at 3 MeV below /7/2. The parameters 
found for the smaller space {N,S,A)= (49,185,0.2), be- 
come (65,285,0.8) for the larger one. The densities coin- 
cide nicely in the region of interest. Discrepancies become 
appreciable after 25 MeV, which provides a preliminary 
indication: binomial densities can be trusted in an inter- 
val of « 0.15S' above the ground state. It is very much in 
the logic of the construction, that knowledge of the level 
density up to a given energy, can be extended to higher 
energies. 
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FIG. 2. Experimental logarithmic level densities for ®'^Ni, 
compared with binomial ones calculated in two spaces 



The first part of the program is now successfully com- 
pleted. 

To understand in what sense binomial and Hamilto- 
nian moments coincide we derive the general form for 
both. The following notation will prove very useful (a is 
an arbitrary array): 



[a]'^^ XI °"i 



ai ^ aj 



(6) 



If ai can take L values, there are L{L—1) ■ ■ -{L — l + l) = 
i*^'^ terms (Boole's notation). If the indeces are ordered, 
Qi < q;2 • • •, we define analogously [a]'^ = [a]li/Z!. 

For the binomial case we start from the standard result 



AKO^J:^^--"- e 



N 



,t{7i~Np) 



(pe«*+(?e-P*) 



pt\N 



K 



- kaIk] 



Ma'=(^) M{t)\t^a=M 



\t=o 



(7) 



To evaluate A^_r-, consider the multinomial expansion 
(Leibnitz rule) for the K-th t-derivative of a product of 
A'' factors ha t 



N 



Y[ho 



[K] 



E^' 



u[ki] . [fc2] 

'h,t 'h,t 
fci! fe! 



(8) 



E-ny- 



E^'n 



[h^^ 



s,h 



{siy^isV 



Y, s^« = K. 



In the second equality, the Ig factors with common k^ = s 

have been regrouped in [/i''*'] ° and then the original or- 
dering among the factors has been relaxed. Notations 
from Eq. (Q). Note also that (obviously) ^^ Is — N. 
Identifying h^^t 



pe"?* + qe p* , we have 



Kt =pm 



s-l qt 



+ (-)V 



-ie-P*l 



(9) 



Calling K = /laLo (^0 = l,/ii = 0), and N - !„ ^ I, 
the moments follow from Eq. (||), and noting that N'--'^ = 
N\/{N ~ ly. = N\/lo\, (the TV! coming from the equality 
of all the ha^t factors), the result is: 






{siy^isV 



s>2. 



The dominant terms maximize /, subject to the condi- 
tions J^s^s = K, J2^s = I < N. Keeping contribu- 
tions up to 0{l/N), the normalized moments {A4k:/(J^ , 
a^ — Nh2) for K = 2k and K = 2k + 1 respectively, are 



Ar(fc) 
N'' 



(2fc-l)!! 



(2fc + l)! 



iN-k + 1) 



6hl 



k^mq 

~9hf 



kh-i 



3VNh: 



3/2' 



(11) 



They will remain dominant as long as fc < v iV. 

To calculate moments of a Hamiltonian, we write 
Ti^ = Cl2a^a,t)^ ■ In the single particle case ha,t = 
maSa, SO the summands commute, and the multino- 
mial expansion is exactly the same as in Eq. (g), with 
derivatives, [s], replaced by powers s, and we have now 
[/i^] 4 = [e''m'^]j. Since m^ — rria, and there are I factors, 
the number operators contribute as Y[i=i i "^a; ' whose 
trace is tvM^ / D'^^^ (it must be proportional to to*^'-*, and 
unity at the closed shell in which the D orbits are full). 
Defining e^ = e^^D~^ the moments for a one body (r = 1) 
Hamiltonian follow as 



D' 



Mr-Kij:^i^^^Y[ 



S,ls 



D(i) 11 (s!)'.^,! 



s> 1. 



(12) 



In the dilute limit (m < D), D^/D'^^^ « 1, and [e"]^ 



;, — 



the notation in MF, whose result 



Eq. (7)] is identical to Eq. ( |12[ ) (the s = 1 terms cancel 
since TC is taken to be traceless). The complete analogy 
with Eq. (|lO|) is obvious. 

To understand how to dispose of the condition m ^ D, 
let us calculate ct^(to), (use in = D — m) 



C£^l 



E^l 



nia + y ]e 


^eprriamfj)^ 


a^l3 


m m(^) 


mm Y" 


D Z?(2) 


- D-l^ 



(13) 



where we have added and subtracted the terms necessary 
to eliminate the restriction a ^ (3 and remembered that 
^^ Ea = 0. Similarly, for the third moment we find 



M 



r=l 



mm{m — m) 

[D-iy^) 



E 



D 



(14) 



Under particle- hole transformation (i. e. m — > to), 
the even moments are symmetric, and the odd ones anti- 
symmetric, which simplifies enormously the calculations. 
E. g., Eq. (|lj) follows from the argument that a'^im) 
must be a two body operator that reduces to the correct 
value for m,,m = 1. More generally, the leading term for 



(10) A^2/c = -^2fc V('''^(''^))'^ must have the form 






-{2k-iy.\ = T{m,k){2k-iy.\, (15) 



which follows from demanding rank 2k, symmetry in par- 
ticles and holes, vanishing for m,rh < k and corrrect 
value at m.,fh — k (i. e. djT^ — kl/D'^^^). The important 
point is that !F{m, fc) « 1 as long as fc <C to, i. e., as long 
as this term remains dominant the moments are those of 
a Gaussian (CLT). When the leading odd moments and 
subdominant even ones are included we end up with an 
expression identical to Eq. ( |ll| ) once N and p have been 
extracted as explained at the beginning. Hence, we have 
proved that the level density for a one body Hamiltonian 
has binomial moments for K < y/N. A strong version of 
the central limit theorem (CLT is a statement about the 
even-iiT dominant term only). 

The near identity between Eqs. (|^) and (||) is an excel- 
lent reason to expect that for single particle Hamiltoni- 
ans the binomial behaviour extends to the full spectrum, 
since Eq. (Q) is a rigorous mathematical result that ap- 
plies in this case. 

For a Hamiltonian of higher rank, H = '^ W^yZ^Zy, 
where Z^ , Zy create and annihilate r particles, the op- 
erators do not commute. Nevertheless, Eq. (O) is of 
use in giving the correct counting: there are (2fc — 1)!! 
ways of contracting Ti^*"' in fc pairs, fc(2fc -I- l)!!/3 ways 
of contracting 7^^'''+^ in fc — 1 pairs and one triple, etc. 
However, the contributions of each term are different. 
In other words: for the dominant term, say, the factor 
J-(rn, fc) becomes extremely complicated. The problem 



was solved by Mon and French |^. Here we give an idea 
of their result. For simplicity assume that 7i is a two- 
body operator {r—2), and stay in the dilute limit. Then 
the variance a'^{m) = {'}i?)m/dm must be 



a^{m) 



7,(2) 



-a\2), a\2) 






(16) 



Now consider (W*), leading to three possible contrac- 
tions, written in MF as nUnU = A ABB + ABBA + 
ABAB. The first two give m(2)cr2(2))V4 = {a'^{m)f, 
but for the last one we have TO'-^-'(cr2(2))2/4, which van- 
ishes at m = 2. The general result is that for m = 2 and 
K = 2k, only tk = ik-i)/^ terms survive, the Catalan 
numbers, i. e., the normalized moments of the semicircu- 
lar density for Wigner's GOE ensemble. As m increases 
the number of surviving terms increases rapidly so as to 
have again !F{m, k) ~ 1. 

It is straightforward to apply similar arguments to the 
other dominant and subdominant terms. Thus, we can 
combine the Mon French analysis and the advantages of 
the binomial geometry, to obtain a strong form of the 
CLT, seen to apply generally to higher rank Hamiltoni- 
ans. 

So far the good news. There are no bad news, but 
a hard question: Why stop the analysis at the low mo- 
ments? As mentioned, it is practically certain that for 
rank one the proof of binomial behaviour must extend to 
the full spectrum. The strong formal analogy between 
rank one and higher should encourage the generaliza- 
tion. The catch is that, more often than not, systems 
undergo phase transitions. As there is little risk in at- 
tributing them to some form of collectivity, we can guess 
that the binomial forms will be valid in the absence of 
strong enough collectivity. Random Hamiltonians fulfill 
this condition, and I propose to include them in a larger 
class: A Hamiltonian acts as random if it does not have 
strong enough collective components. 

There may be two reasons for the good performance 
of binomials in the nuclear case. One is the strong dom- 
inance of the single particle field. The other is the lack 
of sufficient collectivity in the nuclear Hamiltonian. It 
is certainly not random, as it contains sizeable pairing 
and quadrupole forces [|l2[. However, they do not seem 
to be strong enough. A quantitative estimate of the rela- 
tive strength of the different components is given by a^ : 
we have already seen that the single particle contribu- 
tion is far stronger that the two body part, and from |l3] 



we know that in the latter, pairing plus quadrupole only 
contribute a fraction of the total. Sufficient to give them 
a capital spectroscopic role but, apparently not the pos- 
sibility of distorting the binomial forms. 

The formalism is ready to examine the problem, which 
will become even more interesting with the suggestion 
that Hamiltonian matrices at fixed quantum numbers al- 
ways have binomial level densities |l3] . 
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